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■ Based on the solution of a boundary problem for disconnected (Killing) horizons 
and the resulting violation of characteristic black hole properties, we present a non- 
existence proof for equilibrium configurations consisting of two aligned rotating black 
holes. Our discussion is principally aimed at developing the ideas of the proof and 
summarizing the results of two preceding papers (Neugebauer and Hennig, 2009 
[29^1, Hennig and Neugebauer, 2011 (l2l]). Prom a mathematical point of view, this 

Q , paper is a further example (Meinel et al., 2008 0) for the apphcatfon of the inverse 

("scattering") method to a non-linear elliptic differential equation. 

U • 

^ ■ I. INTRODUCTION 

o _ 

The investigation of gravitational interactions in static two-body systems dates back to 
10 ■ the early days of General Relativity and was initiated by Hermann Weyl and Rudolf Bach. In 
i/^ \ a joint paper [2:] they discussed, as a characteristic example, an axisymmetric configuration 
O ■ consisting of two "sphere-like" bodies at rest. Bach, who constructed a corresponding solu- 
tion for the vacuum region outside the bodies by superposition of two exterior Schwarzschild 
solutions, noted that this solution becomes singular on the portion of the symmetry axis 
between the two bodies as expected. In a supplement to Bach's contribution, Weyl fo- 
cused on the interpretation of this type of singularity and used stress components of the 

■ energy-momentum tensor to define a non- gravitational repulsion between the bodies which 
compensates the gravitational attraction. Weyl's result is based on some artificial assump- 
tions but implies an interesting question: Are there repulsive effects of gravitational origin 
which could counterbalance the omnipresent mass attraction? 

Newtonian approximations tell us that the interaction of the angular momenta of rotating 
bodies ("spin-spin interaction") could indeed generate repulsive effects. This is a good 
motivation to study, in a rigorous way, the equilibrium between two (aligned) rotating black 
holes with parallel (or anti-parallel) spins as a characteristic example for a stationary two- 



body problem. In preceding papers [12|, |29| , which involved degenerate ( "extreme" ) black 
holes we came to a negative conclusion. This paper is meant to summarize the steps of this 
non-existence proof, to point out the main points of the matter and to refer non-specialists 
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FIG. 1: Illustration of a two-black-hole configuration with one degenerate (point-like) horizon 
^(1) {g = 0X = Ki) and one sub-extremal horizon Ti^"^^ {q = 0,Ks > C ^ ^i) i^i Weyl-Lewis- 
Papapetrou coordinates. A'^, and A~ denote the three parts of the axis of symmetry. 



to papers dealing with the technical details. It should be noted that a non-existence proof 
for special symmetric equilibrium configurations by Beig et al. jsl, 0] is essentially based 
on symmetry arguments and does not apply to black holes with different horizon areas and 
angular momenta. 

Another aspect with some relevance for the two-black-hole configurations in question is 



the interpretation of the so-called double-Kerr-NUT solution [17|, |2JJ, a seven parameter 
solution constructed by a two-fold Backlund transformation of Minkowski space. Since a 
single Backlund transformation generates the Kerr-NUT solution that contains, by a special 
choice of its three parameters, the stationary black hole solution (Kerr solution) and since 
Backlund transformations act as a non-linear superposition principle, the double-Kerr-NUT 
solution was considered to be a good candidate for the solution of the two-horizon problem in 
question and extensively discussed in the hterature I17H21|. l30l. l33j. However, there 

was no argument that this particular solution is the only candidate. Therefore, defects of 
this special solution would not a priori imply a general non-existence proof for our stationary 
two-black-hole problem. 

In papers we could remove this objection and show that the discussion of 

a boundary value problem for the Ernst equation which represents a part of the vacuum 
Einstein equations, necessarily leads to a subclass of the double-Kerr-NUT solution. This 
result is in line with a theorem of Varzugin (sil . 32 1 which says that the 2A^-soliton solution 
by Belinski and Zakharov contains all possible solutions (if any exist) corresponding to an 
equilibrium configuration of black holes. 

The subclass is characterized by a set of restrictions for the parameters of the general 
double-Kerr-NUT solution. These restrictions ensure the "correct" behavior of the double- 
Kerr-NUT solution along the axis of symmetry and the horizons. For solving the restrictions 
we could go back to the already mentioned discussions of the equilibrium conditions for the 
double-Kerr-NUT solution. After a too restrictive ansatz in [VA , Tomimatsu and Kihara 15 



30| derived and discussed a complete set of equilibrium conditions on the axis of symmetry. 
Reformulations and numerical studies by Hoenselaers jl3] made plausible that the double- 
Kerr-NUT solution cannot describe a two-black-hole equilibrium if the individual Komar 
masses of its two sources ("horizons") are assumed to be positive. Manko et al. [20| and 
finally Manko and Ruiz ||21|] were able to prove the conjecture. The critical point of this proof 
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in view of a non-existence theorem is, however, the presumed positiveness of the two Komar 
masses. To the best of our knowledge there is no argument in favor of this assumption (on 
the contrary, Ansorg and Petroff [1] gave convincing counterexamples). Instead, we replace 
the Komar mass inequality Mj > for each black hole (2 = 1,2) b y in equalities connecting 
angular momenta Jj and surface areas Ai (87r|Jj| < Ai, i = 1,2) These relations are 
based on the causal structure of trapped surfaces in the interior vicinity of any event horizon 
[gI. In the case of two non- degenerate black holes it turns out that one of the two inequalities 
is always violated, i.e. one of the sources cannot be a black hole. 

This type of non-existence proof avoids more laborious investigations of the domain off 
the axis of symmetry. In some degenerate cases, see Fig. [H we need additional eliminating 
criteria such as the positiveness of the total (ADM) mass or the absence of singular rings 
for proving the non-existence. We summarize the results of the discussions of all subcases 
in Sec. ED 



II. A BOUNDARY PROBLEM FOR DISCONNECTED HORIZONS 

The exterior vacuum gravitational field of axially symmetric and stationary gravitational 
sources can be described in cylindrical Weyl-Lewis-Papapetrou coordinates {g,(,'^,tY, in 
which the line element takes the form 

ds^ = e-2^ [e'^'idg^ + dC') + qM^p'] - e^"" (dt + a dipf, (1) 

where the "Newtonian" gravitational potential U, the gravitomagnetic potential a and the 
"superpotential" k are functions of g and ( alone. At large distances r = | a/^^ + I ^ oo 
from isolated sources located around the origin of the coordinate system, r = 0, the spacetime 
has to be Minkowskian, 

r^oo: ds^ = dg^ + d(^ + gMy:>^ - dt^. (2) 

According to the objective of this paper (see Fig. [1]) we will exclusively discuss gravitational 
fields ([T]) under condition ([2]). 

Metric ([1]) admits an Abelian group of motions G2 with the generators (Killing vectors) 

= 61 (stationarity) 
?7* = (5^ (axisymmetry) 

where the Kronecker symbols Si, 5^ indicate that ^* has only a t-component whereas rj^ points 
in the azimuthal (if) direction, r/* has closed compact trajectories about the axis of symmetry 
and is therefore space-like off the axis (and the horizons). ^* is time-like sufficiently far from 
the black holes but can become space-hke inside ergoregions. Obviously, 

e2^ = -fe., a = -e-2V6 (4) 
is a coordinate-free representation of the gravitational potentials U and a. 



^ In the following, we also use the complex coordinates z — g + i( and z — g — t is the time coordinate. 
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According to Carter's theorems [8|] we can assume that the event horizons of the two 
black holes under discussion are Killing horizons. Here a Killing horizon can be defined by 
a linear combination ^' of the Killing vectors ^ and t], 

e = ^+^v (5) 

with the norm 

where i? is the constant angular velocity of the horizon. A connected component of the set 
of points with e^^ = 0, which is a null hypersurface, (de^^, de^^) = 0, is called a Killing 
horizon T-L{C)y 

no-. e2^ = -(r,O = 0, (de2^de2^) = 0. (7) 

Since the Lie derivative C^i of e^^ vanishes, we have (i^', de^^) = 0. Being null vectors on 
^(^')' ^' de^^ are proportional to each other, 

no : de^^ = -2<'. (8) 

Using the field equations one can show that the surface gravity k is a constant on 'H(^'). 

In the Q-C, plane {t = constant, = constant) of the Weyl-Lewis-Papapetrou coordinate 
system ([T]) horizons are located on the C-axis {g = 0) and cover a finite portion of the axis 
(?/(^) in Fig. d]) or shrink to a single point ("H^^^ in Fig. [Tjjsj. It turns out that extended 
horizons ("sub-extremal horizons") and point-like horizons ("degenerate horizons") require 
different considerations. Note that a Killing horizon is always a two-surface in the time slice 
t = constant. The degeneracy to a line or a point is a peculiarity of the special coordinate 
system. 

The dashed line in Fig. [1] sketches the boundaries of the vacuum region: 
are the regular parts of the (^-axis (axis of symmetry), "H*-^-* and "H^^-* denote the two Killing 
horizons (Fig. [T]shows a point-like and an extended horizon), and C stands for spatial infinity. 

The gravitational fields a, k, U have to satisfy the following boundary conditions 



C 



a = 0, k = 0, (9) 
l + i7ia = 0, i = l,2, (10) 
a^O, k^O, (11) 



where i7i and i72 are the angular velocities of the two horizons. Equations (Q characterize 
the axis of symmetry (rotation axis). The first relation originates from the second equation 
in (jlj), since the compact trajectories of rj with the standard periodicity 2tt become infinites- 
imal circles with the consequence 77 — ?■ 0. The second relation is a necessary condition for 
elementary fiatness (Lorentzian geometry in the vicinity of the rotation axis). Equation f IlOp 
is a reformulation of Eqs. ([7]) (e^^ = 0) and since the horizons are located on the (^-axis 
{g = 0); see Fig. [TJ Finally, Eq. ( ITT]) ensures the asymptotic fiatness of the metric ([T]); see 

(ED. 



In our discussion we will essentially use the Ernst formulation of the field equations [10 
For this purpose, we introduce the complex Ernst potential 

f = e'^ + ife, (12) 

where the (real) twist potential b is defined by 

a,i? = ^e"^^6o a,C = (13) 
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In this formulation, a part of the Einstein vacuum equations is equivalent to the complex 

Ernst equation 

m){f,,,+ia+-/o) = fl+fi- (14) 

Obviously, the imaginary part of the Ernst equation is nothing but the integrability condition 
. On the other hand, the condition b ^c = b^g leads back to the field equation for 



^C^?• 



a as an element of the original Einstein equations. 

The metric potential k can be calculated from / via a line integral. 



k. 



Q 



k, 



2q 



U,Uc + -re 



(15) 
implies 



The result of the integration does not depend on the path of integration, since 

Equations f lT^ and f|T5|) are completely equivalent to the Einstein vacuum equations. 
Thus one can first integrate the Ernst equation to obtain e^^ and b or, alternatively, e'^^ and 
a (see ( !T3|) ). and determine the remaining metric potential k (see ([T])) by line integration 
afterward. Taking advantage of this circumstance we will first analyze the boundary problem 



C 



a = 0, 
1 + f2ia 
[/ ^ 0, 



= 0, i 
a 0, 



1,2, 



(16) 
(17) 
(18) 



for the Ernst equation ( fT4l) . Note that the connection between b = '^f and a is non-local; 
see ([131). 



III. THE INVERSE METHOD 



The inverse (scattering) method, first applied for solving initial (value) problems of special 
classes of non-linear partial differential equations in many areas of physics (such as Korteweg- 
de Vries equation in hydrodynamics, non-linear Schrodinger equation in non-linear optics 
etc.), is based on the existence of a linear problem (LP), whose integrability condition is 
equivalent to the non-linear differential equation. Surprisingly, the Ernst equation has an LP. 
This fact is the background for the rich gain of exact solutions with interesting mathematical 
properties; see js], Hif, [13, [22[ [23, [28[ and the references therein. However, we cannot expect 
a priori one of these solutions to solve our physical question, and are therefore referred to 
methods applicable to boundary (value) problems. In a sense, we can try to borrow ideas 
from the above mentioned analysis of initial (value) problems. 



We use the LP 





M 



+ A 



M \ 1 
iV j + A 



N 

M 

M 

N 



(19) 



where the pseudopotential ^{z, 2;, A) is a 2 x 2 matrix depending on the spectral parameter 



A 



iz 



K 



K + iz'' 



(20) 
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as well as on the complex coordinates 

z = g + K, z = Q-K, (21) 

whereas M, N and the complex conjugate quantities M, N are functions of z, z (or g, () 
alone and do not depend on the constant parameter K. Since the integrability conditions 
^^zz = ^,zz niust hold identical in K (or A) they yield the first order equations 

M^z = M{N ~ M) - ^{M + N), N-, = N{M -N)-^^{N + M) (22) 



f f 
M = = (23) 



with the "first integrals" 



where / is any complex function of z, z. Resubstituting M and in Eqs. fl22|) one obtains 
the Ernst equation f|T^ for f{z, z). Thus the Ernst equation is the integrability condition of 
the LP f|T9|) . Vice versa, the matrix ^ calculated from M, does not depend on the path 
of integration if / is a solution of the Ernst equation. 

Without loss of generality the matrix ^ may be assumed to have the form 

^^fi^ieX,^) ^(^,C,-A) \ ^24) 
V x{qXA) -x(^,C,-A) ) ' ^ ' 

Note that both columns are independent solutions of the LP. The particular form of fl2^ is 
equivalent to 

*(-A)=(j_°i)*(A)(';j). (25) 

Furthermore 

V'(f?,C,i) =x(^,C,A) (26) 

due to the special structure of the coefficient matrices of the LP. 
For K ^ oo and A — t- — 1 the functions if), x may be normalized by 

^(f?,C,-l) = x(f?,C,-l) = l- (27) 
For K oc and A — ?■ 1 one finds 

f{Q. C) = x{q. C, 1) {Kq. C) = i^{Q. C, 1)) (28) 

as a consequence of the LP (1191) . Thus one obtains the solution / of the Ernst equation from 
the pseudopotential ^ (x = ^21, ^ = ^11) for a particular choice of the spectral parameter 
A. Similarly, the gravitomagnetic potential a can be rediscovered at — )■ 00 and A = 1 in 
the first derivatives of the pseudopotential # (27| . 

a(f.,C) = -f^e-2^ (|^[x(-A)-^(-A)]) -C. (29) 

K ^ 00 
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To prove this one has to make use of the LP (fT9|) and Eqs. (123|) and (fT3i) . The arbitrary 
real constant C may be fixed by setting a = for r = \/ + C"^ — ?■ oo. An alternative form 
of (El is 



a(f^,C) = ie-^M^'^[x(-A)-^(-A)] 



- C. (30) 

A = 1 ^ ^ 



— >■ oo 



(Note that 8/ OK = X^xd/dX with A from ( |20|) .) The idea of the inverse (scattering) method 
is to discuss ^, for fixed but arbitrary values oi g, ( { holomorphic function of 

A. We will show that the boundary values (19|)-( [TT]) together with the LP (fT9|) yield the 
necessary information that enables us to construct ^{g,(,^)- (The Ernst potential and 
the gravitomagnetic potential can then be determined in a simple way; see ( 128|) . (129|) ). To 
realize the program we will integrate the LP along the dashed line in Fig. [1] starting from 
and returning to any point g = 0, ( G A'^. Obviously, A degenerates at ^ = to A = ±1, 
for which reason we perform the discussion in the two sheets of the Riemann if-surface 
connected with A according to (1201) . Note that the mapping of the Riemann surface of K 
onto the A-plane depends on the parameters g, Thus one has movable branch points 
Kb = iz, Kb = —iz and the branch cut between them changes with the coordinates. 

IV. SOLUTION OF THE BOUNDARY PROBLEM 

A. Integration along the boundary 

Integrating the LP ([19]) along 'H(^), "H^^)^ ^nd using ([23]) one finds for the 
values of ^ on each interval 1,1 = A^,'H^^\A'^,'H^'^\A~, the representation 



f 1 ^ rl rl_( A'{K) B\K)\ 



where is the value of the Ernst potential in the interval I and A\K), B\K), C\K), 
D\K) are integration "constants" depending on K alone. On the axis of symmetry (/ = 
A^, A^) / is a function of C alone {g = 0, ( e A^, A^), 

/ = ^±,^°: f = f\0. (32) 

The same holds on extended horizons which can be characterized hj g = 0, too. Weyl- 
Lewis-Papapetrou coordinates are, however, inappropriate for integrating along degenerate 
■'point-like" ) horizons. As was shown by Meinel, see 2^ , this defect can be repaired by the 



introduction of suitable (local) coordinates in the vicinity of such a horizon. To demonstrate 
the procedure we consider the point-like horizon in Fig. [T] at ^ = i^i and replace g and ( 
by polar coordinates {R, 9), g = Rsin 9 , ( = Ki + Rcos9, in which the horizon is described 
by i? — 0, 6' G [0,7r], i.e. by a line in an R-9 diagram. Performing the integration along 
the degenerate horizon in these coordinates one indeed obtains the structure fl3T]) with 
= f\9). Thus we have 

/ = ^(1) -^(2) . fi ^ |/^(C) for extended horizons ^^^^ 
' rr/r,\ point-like horizons 
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The definition of our Killing horizons is based on the Killing vectors ^' = ^ + ilrj, see ([5]), ([6]), 
([7]), where Q = Q\ for "HS^^ and Q = Q2 for 'HS^^ . To exploit the characteristic properties 
of the horizons such as ([7j) we discuss metric potentials , a! constructed from rj' = r) 
according to (jl]). Thus we obtain the transformations 

(1 - f2a')e^^' = (1 + [2a)e^^, (35) 

where i7 = for n^^^ and ^2 = ^22 for H^^). It can easily be verified that the line element 
retains its form ([1]) after the coordinate transformation 

g = g, C = C, V t' = (36) 

where i7 = Qi,Q2i if one simultaneously replaces e^^ and a by e^^' and a'. Because of 
f p6|) we may call the primed quantities e^^ , a', etc. "corotating potentials" . Applying f lT3|) 
to the primed potentials a', 6' and using flM|) . f l35|) . f|T3|) one obtains the corotating Ernst 
potential /' from / and finally the corotating quantities M', A^' via fl23|) . This procedure 
ensures from the outset that /' satisfies the Ernst equation and guarantees the existence of 
an LP ( fT9l) in the corotating system. The ^-matrices of the two systems of reference are 
connected by the relation 

= T^*, (37) 

where 

^ ( l + Qa-Qge-^"" \ ^ ■ ( ^ ■ \n ( (no^ 

^^=[ l + Qa + Qge-^^)^'^^^''^^'' U ij ^ ^ 

with i7 = This can be checked up by a straightforward calculation. 

The matrix ^ as defined in ^IM can be considered as a unique function of A, which 
is therefore defined on both sheets of the i^'-surface. From this point of view, Eqs. f l3ip 
determine ^ on one sheet only, say, on the upper sheet with A = 1.^ Its values on the 
other (lower) sheet with A = — 1 result from (!25|) . Consider now ^ along C. It does not 
depend on g, (, i.e., *[C] = ^{K), since M and vanish (see ([23]), (^2)). Along C, the 
spectral parameter A ( 120|) changes from A = lat^ = 0,C'~^ootoA = — latf) = 0, 
( — )■ — oo, i.e., starting in the upper sheet at the points of intersection A'^/C, one arrives in 
the lower sheet at A~/C. (Using polar coordinates ( = rcosO, g = rsinO, 9 G [0,7?], one 
has A = ^JJl^^^V^^y[{K'^^^r^^ — > A = e'^ as r — )■ oo and therefore A = 1 for 6* = and 
A = —1 for 6' = TT.) Hence, to return to the upper sheet (A = 1), one has to apply f l2^ to 
^[C] = ^[A'^/C] = lim^_^oo ^[^^l? where we have used that # has to be continuous at the 
point of intersection A'^ /C. Thus we obtain 

*|CM-1 = Ito #1^-1 = ( J ) { Ito ( ; J ) , (39) 

and together with fl3Tl) 

fA-{K)B-{K)\_fOl\f A^K) B^K) \ f l\ . . 

\C-{K) D-{K) )~ \ l Q ) \ C+{K) D+{K) J \ l J ' 



^ A = 1 is arbitrarily ascribed to the upper sheet. 
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where A~ means for I = A~ , etc. 

In a similar way one can use continuity arguments to interlink the ABC D-matrices at 
the other points of intersection. Since the Ernst equation has to hold at all these points, / 
must be continuous and unique there. Note that the intersection values 

/i = /[^+/^«], /2 = /[H«M% f3 = f[AV'H^% U = f[n^^\A-] (41) 

are purely imaginary: As a metric coefficient, e^^ = —g'^^ is continuous at the points of 
intersection. According to flMl) and (Q e^^ = e"^^ on the regular parts of the ^-axis. Since 
e^^ = on the horizons, e^^ has to vanish at these points. 

Furthermore, the validity of the Ernst equations of the non- rotating and corotating system 
at the points of intersection implies that ^ and ^' must be continuous there as well. By 
way of example let us consider the point of intersection A~ /H^"^^ (f? = 0, C = -^4) of Fig. [1], 
i.e., the transition from a regular axis {A~) to an extended horizon ("H^^-*). According to 
(I3T|) and dH]) one has 

^[A-/H^'^] = ( _\ )l-=(^ _\ ) L^'\ (42) 

where L' = for / = A^ and -L*^^-* = for / = 'HP'\ Note that the determinant of the 
first factor of the matrix products vanishes such that one may cancel one line of the matrix 
equation. The corresponding equation for at 2; = iK/^ follows from Eqs. (!3T|) . (l37Il and 
dMl) under conditions ([9]) (a = on A-) and ^ (1 + i^aa = on 

Again, the matrix equation consists of two identical lines. Combining a line of fH2|) with a 
line of (H3!l one obtains 

( /4 + 2iQilK - 7^4) -1 ) ^ = ( 2iQ2{K- ^4) ) 

Consider now two extended horizons: H^^^ : C e [Ki,i^2], "H^^^ : C e [i^3,K4], Ki> K2> 
> K4. Following the idea that led to equation (jUj) one can continue the connection of 
adjoining L-{ABCD-) matrices. Involving (HO!) and (jHj) and defining 

Zp /J' i = 1,2,3,4, (45) 

one arrives at the following chain of equations that refiects the integration of the LP along 
the closed contour ^+H(iU°H(2)^-C^+, 



2iQ^^){K - Ki)) ' V 2il2(2)(ir-K2 

^ ^2in(^){K-K^)J ' V 2if2(4)(ir-ir4)' ' ^ 

= ' 1 ) ( 1 
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where Q^^^ = f?^^) = f2i, f?^^) = ^(4) = 

Point-like horizons can be involved without any difficulty by setting Ki = K2 or/and 
= K4 in relations ( l46l) . Consider, for example, the point-like horizon in Fig. [H Though 

the horizon is placed on a point (as a peculiarity of the Weyl-Lewis-Papapetrou coordinate 

system), the Ernst potential has different values at /Ti^^^ and Ti^^yA^: fi 7^ f2', see (HTj) . 

Setting ( = Ki = K2 does not require special considerations such that one can adopt the 

discussion of the "extended" case step by step. 

Eliminating L^^\ L^, L^^\ by means of (H6l) and defining 



i=l 



where Qi = n^^^ = Q2 = ^^^^ = one arrives at the final result of the integration 
along the boundaries 

(48) 



1 

1 



which specifies the holomorphic structure of the "integration constants" L as functions of 
the complex spectral parameter K. Since the trace of the left hand side of (jUj) vanishes, 
has to be trace free, 

tril+ = 0. (49) 

As a condition identical in K, equation (H^ yields four constraints among i72; Ki — K2, 
K2 - Ks, K3- Ki, fi, U By way of example, 

as a consequence of limx-j^oo tr(i?''^K^) = 0. The asymptotic behavior of is prescribed 
by dSD, (EZD, ([2SD and dSH), where one has to choose I = A+, 



Nevertheless, equation (1481) is not sufficient for determining uniquely. To illustrate the 
degree of freedom, we factorize L^, 

+ _ / F(K) \ (a{K) P{K)\ 
- \G{K) 1 ) \ P{K) a{K) ) ^^^^ 

which is always possible for det Iy+ 7^ and 7^ ±5+; see (IHTj) . Inserting f l52p into 
one obtains 

-G F 
G 



-- R+, (53) 

i.e., F{K) = F{K) and G{K) = —G{K) are uniquely determined. Both functions are 
regular everywhere in the complex i^-plane with the exception of simple poles (two extended 
horizons) or /and confluent poles of second order at most (one or two point-like horizons). 
In the next section, we shall determine the axis values of the Ernst potential from 

and show that a and /3 do not affect these values, i.e., ^ ^ ^ j is a gauge matrix. 
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B. Axis values of the Ernst potential 

The elements of the matrix ^ = dM]) are unique functions of A in the vacuum region 
outside the horizons. That imphes that its elements ip and x must be unique (^n = <Pi2, 
^21 = —^22) at the confluent branch points^ Kb = Kb = ( of the Riemann i^-surfaces 
belonging to axis values g = 0, ( & A'^,A'^, i.e., according to ( l3Ti) . ( 132|) 



* -- 
Then 



/ = ^±,^0; K = Kb = C: 

^ i^\_(rp i\fm)B\o\ ^(^oi\ 0^ (54) 



^ • /K) + A^(C)-s^(C)' ^ ^ 

i.e. one can express the values of the Ernst potential on the regular portions of the axis by 
the "integration constants" A, B, C, D. According to (1521) . this means for f\C) = /^(C); 

or ^ _^ 

^^"-^^ fHO + fHCY ^^^^ ^ fHC) + fHCY ^^^^ 

i.e., the matrix ^ j does not affect the axis values of the Ernst potential on A~^. 

Because of the successive transformation P6l) this holds for /^(C) and /°(C), too. (Note 
that (155!) and ( HOj) imply f'^{C) = 1//^(C), C ^ where and are continuations 

of /+(C), C e A+, /-(C), C e via A(K), B{K), C{K), D{K).) Thus we obtain the axis 
values /^(C), /°(C) via and (15^ from R^{K = (), see ( H7|) . as well-defined functions 
of C. The parameters entering these functions are restricted by the constraints tri?+ = 0. 

Since the matrix ^ j does not affect Z"*^, /° we may set a(C) = 1, /3(C) = with 

the consequence 

[/3{K) a{K))-[0 ij- ^^^^ 

(Note that a{K), /3{K) are analytic continuations of a(C), /9(C)-) The particular choice flSSj) 
is closely connected with a gauge transformation of the matrix ^: Any transformation 

^new «'old ^^^^ ^^^^ ^ lOyj 

with 

a(if) = a{K), h{K) = -b{K); a{K) 1, h{K) ^0 as K ^ 00 (60) 

leaves the normalizations (!24l) - (l27j) unaffected and can be used to remove the a-/3-matrix 
from fl52l) . Hence (159!) is a gauge transformation and one can adjust the gauge so that 



/ F{K) \ / a{K) (5{K) \ 
\G{K) 1)^ [piK) a{K)) 



(61) 



—iz Q as Q ^ Q. 
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in accordance with (l58l) . By (13T!) one has 

J = ^+: ^i2 = V^(-l) = l, <?22 = -X(-I) = 1, (62) 

i.e., the gauge ( l6T|) is equivalent to the formulation of special initial conditions for ip and x 
at some starting point ^ = 0, C = Co ^ -^^5 A = — 1 (if in the lower sheet) of the integration 
along the closed dashed line in Fig. [1] (which we performed with unspecified integration 
constants A+, S+, C+, D+). 

Let us summarize the results of the integration of the LP along the boundary and the 
determination of the axis values /^(C), /°(C)- With the standard gauge ( 1611) and the repre- 
sentation (pTl) we obtain for ^ on 

i-A^: *^(^:^^\).^ ^^-izw)^ 

where F{K), G{K) are elements of the matrix R^; see (1531) and (l47j) . 
must be trace free, 

tril+ = 0, (64) 

see (l49l) . which affects via F{K) and B{K) the constant parameters entering the Ernst 
potential /^(C) on the axis 

I = A^: r(C) = ^^|p, (65) 

see (l56l) . The particular form of the pseudopotential ^ on the intervals 'H^^\ A^, 'H^'^\ A~ 
and the axis values f^{C), /~(C) result from ( l46i) with the "starting matrix" X"*" as chosen 
in ( 16T1) and expressions ( l55l) . 

/^(C) seems to be a quotient of two normalized^ polynomials of fourth degree; see (l56ll . 
(!53|) . (1471) . However, the constraints (149|) take care that the numerator as well as the denom- 
inator are of second degree. Inserting the first equation in fl5^ (/ = A'^) into the second 
one and using (HSj) one obtains 

with the consequence 

[i^^(C)-l](^+^(^)) =0, [i2+(C) + l](_jl(^)) =0. (67) 
By definition (H7I) . the elements i?^ of the matrix obey the conditions 

-^11 ~ -^22 ~ ~^22^ -^12 ~ -^12' -^21 ~ -^21- (^8) 

Hence the two equations fl67|) are complex conjugate. /"""(C) can now be calculated from (167|) 
provided that det(i?"'" — 1) = holds. By definition, det = — 1, such that det(i2"'" — 1) = 
— tri?"*". Hence, 

det(il+ - 1) = 4^ tr = (69) 



^ The fourth order coefficient is equal to one. 
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are alternative formulations of the constraints. With the aid of the polynomial matrix 



1=1 



f~^{C) takes the form 



(70) 



/"'(C) 



1=1 



11 



11 



22 



0. 



(71) 
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According to fl69l) . the constraints r^^ + r = can be reformulated to give 



n(c - 
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1=1 



1=1 



rVnr 



12' 215 



(72) 



where by definition the two factors ( "brackets" ) on the left hand side are normalized complex 
conjugate polynomials of fourth degree {ff^ = —r^i) whereas and r^i are normalized real 
polynomials of fourth degree. 

Identifying the zeros on both sides of fl72|) we see see that each factor on the left hand 
side has to have two zeros of as well as of r^^. (Note that the brackets are complex 
conjugate.) Hence the numerator and the denominator of as factors in fl72|) have two 
common zeros such that has to be a quotient of two polynomials of second degree, 



r(c) 



MO ^ C + qC + r 

4(0 e + sc+f' 



(73) 



where the complex constants q, r, s, t are restricted by constraints (jUj). Axis values of the 
form ([72]) are characteristic for the Ernst potential of the double-Kerr-NUT solution. Before 
continuing f~^{C) to all space, we will discuss the gravitomagnetic potential a on A^, A^, 
-^(i/2)_ Pqj^ ^]-^jg purpose we express x(— A), ■?/'(— A) in f l5U]) by the elements of the last column 
of * in on])- The result 

d 



-2iK' 



:B\K)\ 



c 



(74) 



dK ^"^i^^oo 

tells us that the gravitomagnetic potential has constant values on all intervals A^, A^ 
1-L^^^'^\ Since a^gX) is only determined up to an arbitrary constant, we adjust C so that 



0. 



(75) 



To connect the values of a on adjoining intervals we make use of the successive transfor- 
mations (146 p but forgo for the moment the special gauge (16 ip . Considering the transition 
we assume that the asymptotics of can be described by a power series in l/K, 

with 

D+{K) = l + o{^^^, B+{K) = o(^-^y (76) 

Obviously, our standard gauge (16T1) satisfies this assumption. Applying (1741) to the 1-2 
element of the matrix equation 



1 - 



2i[2W{K -Ki) 



(77) 
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see (H6|) . one obtains 



a(^) = a+ - ^. (78) 



Continuing the procedure according to fj46l) one arrives at 

,o^,a, + _L. a'^)=a»-^. a- = a''^> + ^. (79) 

In our context, i?^^) = f^^^) = j^^^ ^7(3) ^ ^(4) ^ However, the form of Eqs. ([78]), ([79]) 

remains unchanged repeating their derivation with il^^^ ^ f]^^^ ^ f]^^^ in (1^6]) . We 

shall make use of this later on; see (|93]1 . 

We can now make sure that our standard representation (l63l) satisfies the boundary 

conditions ([TB]) and (ITT)) . Indeed, L]^2 = B^{K) = and = ^{^) = 1 obey assumption 

([in]) such that Eqs. ([TS]), ([ZH]) hold. With a+ = (see ([TS])) and f^^^) = i?^^) = j7^, 
^(3) ^ ^(4) ^ ^j^gy ^g^j^g fQ^j^ 

a± = 0, a° = 

a» = -^, . = 1,2 

in accordance with (|T6l) . ([T7|) . 



C. Equilibrium conditions 

We have shown that the Ernst potential ( |7T1) can be reduced to a quotient of two poly- 
nomials of second degree by means of the constraints rn + r22 = 0. However, the explicit 
determination of the coefficients q, r, s, t in ( 175]) turns out to be a subtle point. It depends 
on a suitable reparametrization of the constants i72, Ki — K2, K2 — K^, — K^^^ fi, 
. . . , /4, which permits an easier handling of the constraints. For this reason we introduce 
the functions 

«(C) = ^, aa = l, /3(C) = (3(3 = 1 (81) 

and discuss their behavior at the points ( = K^, {i = 1, . . . ,A) which fix the positions of the 
horizons. To begin with, we examine the configuration of two extended horizons 

n^^^ : Ki>C>K2, n^^^ : K3 > ( > K^, > K2 > K3 > K^. (82) 

Introducing the parameters 

ai = a{Ki), aiai = l, i = l,...,4, 

A = /3(i^.), M = l, ^ = 1,...,4, ^ ^ 



we obtain from ([81]) the two linear algebraic systems of equations 

d2iKi) - aMKi) = 0, n2iKi) - /3MKi) = 0, . = 1,...,4 (84) 

for s, t {s, t); q, r (g, f) with Ki, ai, (3i as coefficients. According to ([FT]) . ([S5]) and ([TTI) and 
using f]*"^ = — r^]^ we have 

g2C7+ ^ {C - K{){C - K2){C - K,){C - K^) ^ ^g^^ 
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where is a real normalized polynomial of fourth degree in C. Frome2^^(-^') = (r+iKi) ^ 
0, 2 = 1, . . . , 4) we get 



with the consequence 



(86) 



(87) 



Hence, can be expressed in terms of (and Ki) alone. Solving the linear equations ( 184|) 
for q, r, s, t and plugging the result into (17^ we arrive at a determinant representation of 
the axis potential on A'^, 



no 



1 Kf 

1 aiK,{C-K,) 



aiiC-Ki] 
1 



Kl 
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Kl 
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Kl 
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K! 



^2 



^3 










aiKi{C-Ki 
Kl 



a2K2iC-K2) a^KsiC-Ks) a^K^{C-K^) 



It can easily be seen that 



/(^^,C) 
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- Ki)^ + g^>0, 


z = l, 


..,4, 





(88) 



(89) 



K-K,li = i. 



(90) 

4 in dHHD by 



is a continuation of f^{C) to all space. (Replace ( — Ki 
r„ i = 1, . . . ,4.) 

f{g, C) is a solution of the Ernst equation [l7|. As we have already mentioned, the other 
gravitational potentials k, a {e^^ = 3?/!) can be calculated from f{g,C) via line integrals. 
This solution of the Einstein vacuum equations is known under the name of double-Kerr- 
NUT solution. Since it can be shown by the inverse scattering methods that the axis values 
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/"•"(C) uniquely determine the Ernst potential everywhere in the g-( plane, f{g,C) as defined 
in (189|) is the only solution of the Ernst equation to the boundary values (IHHl) . Hence, the 
solution of the two-horizon problem must be a (particular) double-Kerr-NUT solution. 

The double-Kerr-NUT solution itself is a particular case (A^ = 2, Minkowski seed) of a 
class of solutions generated by an iV-fold Backlund transformation from an arbitrary seed 
solution of the vacuum field equations 23-25[^. The interrelationship with a class of solitonic 
solutions discovered by Belinski and Zakharov jsl is discussed in 16 



30| . Applying 



Many attempts have been made to establish a connection between the double-Kerr-NUT 
solution and two-black- hole equilibrium configurations [§, Q, 15, 17, 20, 21, 
the boundary conditions 

a = 0, k = , . 



C : 



to /(f?, C)) Tomimatsu and Kihara |15l . |30| derived a complete set of algebraic equilibrium 
conditions on the axis of symmetry connecting the parameters a^, Ki {i = 1,...,4) be- 
tween each other. Particular solutions of the algebraic system involving numerical results 
were discussed by Hoenselaers [13], who came to conjecture that the double-Kerr-NUT solu- 
tion cannot describe the equilibrium between two aligned rotating black holes with positive 
Komar masses. Hoenselaers and Dietz jo], |l3| and Krenzer [l^ were able to prove this con- 
jecture for symmetric configurations Ki — K2 = K3 — K^, fli = The explicit solution of 



21|. Finally, 



the Tomimatsu-Kihara equilibrium conditions was found by Manko et al. [2L ^ 

Manko and Ruiz 2l[ were able to prove Hoenselaers' conjecture. The results derived by 
Manko and collaborators are important steps toward a non-existence proof. In particular, 
we will make use of their solution of the equilibrium conditions. Before using these results, 
we had to formulate and analyze a boundary problem for disconnected horizons, cf. (sij . 
since a non-existence proof cannot be based on an arbitrarily chosen solution, even though 
this solution (here: the double-Kerr-NUT solution) seems to be a promising candidate. For- 
tunately (or, as expected) the analysis of the boundary problem led to the double-Kerr-NUT 
solution. 

There is another critical point in the argumentation of Hoenselaers and Dietz and Manko 
et al. To the best of our knowledge there is no argument in favor of the positiveness of the 
individual Komar masses of interacting bodies and black holes. On the contrary, Ansorg and 
Petroff pLi] have given convincing counterexamples. We replace the Komar mass inequality 
("positivity of the Komar mass of each black hole") by an inequality connecting angular 
momentum and horizon area 
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see Sec. IVII This relation is based on the geometry of 
trapped surfaces in the interior vicinity of the event horizon. 

Our next goal is the reformulation of the constraints (!69|) in terms of the new parameters 
(i = 1, . . . , 4). In order to utilize the tools of the inverse (scattering) method we make 
use of the pseudopotential ^(^, C, A) as a solution of the LP (fT9|) . whose coefficients N{q, (^), 



^ It should be noted that these solutions can also be written as a quotient of two determinants in complete 
analogy to ([Ml); see [l3 or 0. 
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M{g, Q are explicitly known by (189|1 and ( 123|) . The integration of the LP yields (cf. |26| ) 
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(92) 



where = V(C^^^iF + ^ > and i = 1, . . . , 4. Note that the remaining; elements of 
^ can easily be constructed from xioX^^) = '^21- ^11 = ^(^?, C;-^) = xioX^^'^)^^ see 
fl26|) . = C, — A), ^22 = — C) ~'^)- The straightforward verification of fl92|) by 
inserting ^(^), A) into the LP f|T9l) with coefficients fl23|) and fl89|) is laborious. Instead 



26l |. one realizes that ^ z^ , for fixed values of g, ( is regular in the A-plane with the 
exception of a simple pole of first order at A = 00. From Liouville's theorem one may 
conclude that ^ ^^^^ = P + XQ, where the 2x2 matrices P{g, (), Q{g, C) do not depend 
on A. As a consequence of fl2^ . P becomes diagonal and Q off-diagonal with Q12 = Pn and 
Q21 = P22 such that one arrives at the first equation f[T^ . From this equation one obtains 
the second equation flTIJl) by using fl2B]) . Note that the Ernst potential (15^ follows from 
(I92|). /(^, C) = x{q, C, 1) ^ 00) as expected; see (l28|) . 

Having discussed some implications of the reparametrization for the representation of 
the Ernst potential and the pseudopotential ^, we will now formulate the constraints in 
terms of at and Ki, i = 1, . . . ,4. In particular, we will show that representation (!88l) of 
together with the boundary conditions (fT6|) . (fTTl) for the gravitomagnetic potential a 
is equivalent to representation (I7T]) with trr+ = (<(=^ tr = 0). 

Denote the pseudopotential ^ constructed from ^% by and equip all quantities 
derived from with an index a. Since is an integral of the LP fll9p . its values on 
the axis and horizons have the form fl3Tl) . They are continuous at the points of intersection 
A+/n^^\ H(^V^°> ^7^(2), ^(2)/^'. This may directly be verified in ([92]) and for the o ther 
elements of the matrix ^q,. The introduction of rotating systems of reference in IIV Al and 
corotating quantities such as the corotating pseudopotential ^' are closely connected with 
the definition of the Killing horizon, which has a well-defined angular velocity in contrast 
to the intervals [Ki,K2\ and [ir3,i^'4] as "potential" horizons. However, one can introduce 
different corotating systems of reference at the ends of each "potential" horizon with angular 
velocities i7'^*\ 2 = 1, . . . , 4, as defined in fl75]) and fl7^ . 



,(2) 



]'\ ^'''=[a--a^^^Y\ (93) 



^ Due to pp]) . K IS a. rational function of A and can be replaced in all elements of ^ to obtain ^{g, Ci X) as 
a polynomial of fourth degree in A. A — !■ 1/A implies K K and vice versa. 
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where, in general, 7^ 7^ '^a ■ Thus one obtains 
1 



1 



i=l 



ly 



1 



(94) 



where fi = f{g = 0,( = Ki), z = 1, . . . , 4, entering Fi according to f HSj) . can be taken from 
f l89|l . A direct consequence of these equations is 



which may be compared with the constraints 
satisfied if f]^^^ = f]^^\ f^^^) _ ^(4) 







(95) 



49|) . Obviously, the constraints tri2+ = are 



f2o 



= = [a- - aL^)] 



-1 



(96) 



Conversely, we have shown that the representation (!7T|) of f'^{C) together with the con- 
straints trr"*" = = tri2+ implies conditions (IHOl) . Thus we may conclude that the (four) 
conditions (!96l) are a reformulation of the (four) constraints According to (fT6|) . (fT7|) 

the conditions ( !96l) are necessary conditions for the equilibrium of the two-black-hole con- 
figuration. Hence, constraints ( IMll and restrictions (l96l) are equivalent formulations of the 
equilibrium conditions to the respective Ernst potentials. 

The explicit form of the axis equilibrium conditions (l96ll can be evaluated by using (l29l) 
or ( 17^ . In both cases, one may start with expression (1921) that determines all elements of 
(e^^ = 3?/ in fl29|) can be taken from fl89l) ). Equation fl29l) yields the gravitomagnetic 
potential a{g,() everywhere including the axis intervals / whereas equation ([71]) directly 
leads to the axis values of J I = A^, v4°, 'H*-^''^-' {B\K) may be taken from representation 
fl3T]) ). Straightforward calculations result in 
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where 
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with 



n (1 + A^«. 

n=l 

1 + Xiak 



X 



\Kk-(\ 
Kk-C' 



(99) 



^ Henceforth we omit the index a. 
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Note that = ±1 marks the interval /, e.g., I = A^, < K3 < ( < K2 < Ki. 
(AO) = (1,1, -1,-1). 

Inserting determinant expressions (1971) in the axis equihbrium conditions ( p6|) ^ one obtains 
the angular velocities of the horizons l-L^^^"^^ in terms of the parameters ttj, Ki^ i = 

1, . . . ,4 and two conditions (a+ = ao, = a^) that restrict the choice of these parameters. 
These restrictions must be taken into account when examining the Ernst potential (l89il . 

D. Discussion of the solution 

We have shown that the solution of the boundary problem ( fT6i) - (|T8l) for the Ernst equation 
iHM is given by the Ernst potential (l89l) . whose constant parameters ai, Ki, i = 1, . . . ,4 have 
to satisfy the first set of the axis equilibrium conditions ( |96l) . a+ = a", a~ = a° with a^, a° 
from fl97|) . (i?!, can be calculated straightforwardly.) The appearance of the equilibrium 
condition is not promising and makes a comprehensive discussion of the solution difficult. 
We will only list a few aspects of the interpretation. 

1 . Number of parameters 

Written in dimensionless coordinates such as 

= -^, C = ^^, K2s = K2-K3, (100) 

the Ernst potential contains four free parameters. Since the quotient of determinants (!89|) 
remains unchanged under a translation of the two-black-hole configuration along the (- 
axis and a multiplication of g, ( and by a common (real) factor, / depends on the 
coordinates g, ( and the six parameters ai, . . . , 04; K12/K23, K^i/K23, where we have used 
the abbreviation 

Kij = K,-Kj, i,j = l,...,4. (101) 

It can easily be seen that the two conditions = a°, a~ = aP can be rescaled to be written 
in terms of the six parameters alone. Hence / = /(^, Q is a four parameter solution. 

2. Singularities outside the horizons 

The solution is a necessary consequence of the integration of the LP along the boundary 
(closed dashed line in Fig. [T]). We have (implicitly) assumed the validity of the Ernst 
equation everywhere in the enclosed g-C, domain. It is, however, not clear, whether the 
Ernst potential is really free of singularities there. As matters stand at present the question 
must remain undecided. Interestingly, the solution of the static two-horizon problem on the 
level of the Laplace equation AU = (which is the static form of the Ernst equation) has 
no singularities outside the horizons. (Applying the inverse formalism one simply obtains a 
superposition of two Schwarzschild solutions in Weyl-Lewis-Papapetrou coordinates). The 



We omit a as arranged. 
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so-called "conical" singularity of the metric on that forbids the existence of static two- 
black-hole configurations only appears when the metric coefficient e^^ is involved. Calculated 
from a regular solution of the Laplace equation (the "double-Schwarzschild" solution), e^'^ 
violates the regularity condition e^*^ = 1 on A^; see (jH])- We will discuss this condition (which 
ensures elementary flatness on the axis of symmetry) in the next section. 



3. Generalizations 

Obviously, our analysis of the two-horizon problem can easily be extended to an arbitrary 
number n, n > 2, of aligned disconnected horizons. Integrating the LP (I19p along the 2n-|- 1 
intervals / (n horizons H and n + 1 "regular" intervals .4) one arrives at a representation 
of the form (I3T1) for each of the 2n + 1 intervals. Consequently, one has 2n + 1 matrices 

and 2n matrices i^,; cf. fH5|) . Replacing the symbol Y[t=i P71) by Hi^i '^^^ obtains, 
via fl53|) and by fl56|) . a representation for the Ernst potential f~^{C) on A~^. The rational 
structure of this potential (/^ is a quotient of two normalized polynomials of equal degree 
in () is a characteristic feature of solutions to the Einstein equations derived by iterative 
Backlund transformations of the metric of the Minkowski space [l6| , (2^. or, equivalently, 
by the Belinski-Zakharov approach [sl]. This confirms Varzugin's result j3l| which says that 
any equilibrium configuration of aligned black holes can be described by a Belinski-Zakharov 
solution [sj. After a reparametrization in full analogy to flHTl) . /(f?, C) again turns out to be a 
quotient of (2r;,+l) x [2n+l) determinants whose structure is an obvious generalization of the 
determinants in (l89l) . Finally, the equilibrium conditions can be derived straightforwardly. 



V. ELEMENTARY FLATNESS AND EQUILIBRIUM CONDITIONS 

We know from the Bach-Weyl paper [3] that the metric coefficient e^^ is a measure for 
the interaction of the two black holes. To guarantee equilibrium, k has to vanish on the 
portion of the axis of symmetry between the two black holes. From a geometrical point of 
view, the condition e^^ = 1 on ^ is a necessary condition for elementary fiatness (Lorentzian 
geometry) of spacetime in the vicinity of the rotation axis A. Our discussion of the metric 
potential k = k{g,() is based on Kramer's representation [l^, which is a result of the 
integration of the defining relation (HM with f{g, () from (15^ . It turns out that k{g = 0, () 
is a step function with constant values on the intervals A"^, A'^. In particular, one has 
g2fc+ _ g2fc ^ Equation (1151) determines k up to an additive constant, which can be chosen 
such that k^ = with the consequence e^^ = 1. For this choice, e^^" takes the form 

with , I = A'^,A'^ as in fl98|) . Obviously, the equilibrium condition e^'^" = 1 has two 
solution branches, = ±H^. Tomimatsu and Kihara jsof ruled out the condition = 
[H^ = together with a"*" = a~, a~ = leads to overlapping horizons). = —H^ 
yields 

Q!lQ!2 + a3Q!4 = (103) 
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and effects a considerable simplification of the equilibrium conditions (p6ll . Putting ( I97|) -( p9|) 
in = aP , a~ = aP and observing = —H^ one finds 



03(1 - a^fiK^^i + K32) - (1 - a3a4)a43^3i^32 = 
ai(l + 02)^(^41 + ^32) + (1 - aia2)a2i^3i^4i = 0, 



(104) 



where 



Q!,: — Q!j 



ir., = K.-ir„ a., = -^i— ^. (105) 

For point-like horizons Ki = K2 or/and K3 = K4, the parametrization ( l83l) does not ap- 
ply, since the mapping of the four coefficients q, r, s, t in f l75]) onto less than four parameters 
ai (cti = 02 or/and 03 = 04) is not invertible. However, the invertibility can be restored by 
introducing the derivatives a'(C), /3'(C) iii the confluent points Ki = K2 or/and K3 = K4^. 
As was shown in [l2| this concept makes it possible to consider the equilibrium conditions 
for configurations with degenerate horizons as particular cases of (11031) . (I104p : To describe 
point-like configurations one has to set 

K2 = Ki, a2 = ai, ^21 = a'(-K'i) = -iai7i or/and 
K4 = Ks, 04 = as, ^43 = a'iK^) = -1^373, 

where 71 and 73 are real constants. 

To introduce ^21 and 043 in the Ernst potential (189|) . one has to subtract the second 
column from the third one and the fourth column from the fifth one and to decompose the 
a-parameters into symmetric and antisymmetric parts (e.g. ai = |(q;i -|- ^2) + |(ai — 02), 

0:2 = |(ai + «2) — |(«i — a2))- 

We are now prepared to discuss all constellations (extended/extended, extended/point- 
like, point-like/point-like) by using ([89]), f lT03|l and ffTM]) . 



VI. BLACK HOLE INEQUALITIES AND SINGULARITIES 

A. Two sub-extremal black holes 

We start the discussion of the different types of possible two-black-hole equilibrium con- 
figurations by considering spacetimes containing two black holes with extended horizons. 
Following Both and Fairhurst we will assume that a physically reasonable non-degenerate 
black hole should be sub-extremal, i.e. characterized through the existence of trapped sur- 
faces (surfaces with a negative expansion of outgoing null geodesies) in every sufficiently 
small interior neighborhood of the event horizon. As shown in [ll| , the presence of trapped 
surfaces implies the inequality 87r| J| < A between angular momentum J and horizon area 
A of the black hole. In a regular spacetime with two sub-extremal black holes, both black 
holes have to satisfy this inequality individually, 

87r|Ji|</li, i = l,2. (107) 

This is the key ingredient for the non-existence proof of two-black-hole equilibrium config- 
urations, as we will see below. 

In order to test these inequalities for that subclass of the double-Kerr-NUT family of 
solutions which describes two gravitating objects with extended horizons, we have to solve 
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the general axis regularity conditions fllOSp . f ll04p . For Ki > K2 > > these conditions 
can be written as 

aia2 + a^a/i = (108) 

and 



04)2 2 (l-"3)^ . lKiiK2A 



-w = , w:= J € [1,00) 

"4 as V K13K23 



"2 "1 V Ki^Ku 



w'a^ + iea o? + \w'ea 

Oil = — ; : , "2 

w' — lea 



wa^ — a — wa 



as = , ^4 



w — a 1 — wa 

where a := \/—olxOL2 = y/a^al =: e"'^, aa = 1 and e = ±1. 
Now we are in a position to calculate the quantities 



(109) 



As shown by Manko et al. |21;|], these equations can be explicitly solved for ai, . . . , a^, 

(110) 



Pi-=^^ i = l,2, (111) 

in terms of the parameters w, w' and G [0, 2tt). The result is 

1 + <Pw' w(w , ^ 

Pi = £^^ 7, P2 = e— 112 

with 

:= COS0 + £:sin0, e = ±1. (113) 

Rewriting the inequalities in fll07p as pf < 1, i = 1,2, we find with the previous formulae 
the two conditions 



w 
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2$w' + l<0 and w^-2$w + l<0. (114) 



Since the latter inequalities imply $w' < and $w > 0, we arrive at a contradiction, because 
w and w' are positive. In this way, we have shown that two sub-extremal black holes cannot 
be in equilibrium. 



B. One sub-extremal and one degenerate black hole 

Without loss of generality, we may assume that the upper horizon is point-like [Ki = K2) 
and the lower one is extended {K2 > > K4) as sketched in Fig. [H Then, the equilibrium 
conditions fllOSp . f ll04p can be written as 

af + asa/i = 0, (115) 



with the solution 



«! = iea, a^a^ = o? , 



(116) 
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where e = ±1 and a G C, aa = 1, and 

(1 — 04)^ = , (a — ie) [yiK23w{ea + i) + i{w + l){ea — i)] = 0. (117) 

a 

As solution to these equations we obtain two sets of parameters. In the first solution branch, 
the parameters ai, 71, 0:3 and 04 have to be chosen according to 

i(w + l)i — ea wa^ — a — wa /-,-,r.\ 

ai = lea, 71 = — — , a^ = , a^ = , (118) 

WK23 1 + ea w — a 1 — wa 

which is the limit Ki — )■ K2 w' — t- 1) of solution flllUp for extended horizons. This 
family of solutions depends on the two parameters a and w (and on two additional scaling 
parameters, e.g. Ki and -^'23)- 

The second solution branch of the equilibrium conditions is given by 

TT, l-iew 1 + iew 

tti = -1, 7i e R, as = z——. , 04 = -:; : , (119) 

1 + lew 1 — lew 

i.e. the corresponding Ernst potential depends on the two parameters 71 and w (plus two 
scahng parameters). Interestingly, this solution has no counterpart in the case of extended 
horizons. 

The desired non-existence proof follows the same idea for both families: The ADM mass 
M of the spacetime can be expressed in terms of the two parameters (and the additional 
scaling parameter K23). The resulting expressions can be estimated, using the inequality 
87r| J| < A for the sub-extremal object. We obtain 

M < (120) 

in contradiction to the positive mass theorem. This indicates the presence of unphysical 
singularities and we conclude that configurations with one degenerate and one sub-extremal 
black hole cannot be in equilibrium. 



C. Two degenerate black holes 

In the case of possible equilibrium configurations with two degenerate black holes^, we 
find three one-parametric families of candidate solutions: The equilibrium conditions are 
now 

al + al = 0, (121) 

which is solved by 

ai = iea, a^ = —a, (122) 
with e = ±1 and a G C, aa = 1, and 



(a 



+ 1) [{a - 1)73/^23 - 2i(a + 1)] = 0, {a- ie) [7i^23(£a + i) + 2i{ea - i)] = 0. (123) 



See also [7j for a discussion of properties of spacetimes with two degenerate objects. 
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Equation f ll23p has three different solutions. The first one is 

2i 1 + iea 2i a + 1 

ai = lea, 03 = -a, 71 = — ■ — , 73 = — (124) 

A23 1 — lea K23 a — I 

and depends on one free parameter a. This solution can be obtained in the limit Ki — )■ K2, 
K^i^w-^ 1, w' ^ 1) from f lTTU]) . 
The second and third solution branches are 

«! = -ie, as = 1, 7i = J^, 73 G R, (125) 



K. 
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and 



2e 

ai = -1, 03 = -le, 71 e R, 73 = (126) 
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where now 73 or 71 are free parameters. 

In order to show that the above solutions do not lead to regular two-black-hole configu- 
rations, we calculate the ADM mass M. For the first solution branch, the result is 

M = ^. (127) 

3 + y2cos(£0 + f) 

Obviously, M is always negative and we arrive again at a contradiction to the positive mass 
theorem. 

For the second and third solution branches, M has the form 

M = -^-^^, 73 = 5737^23. (128) 
2 73-1 

Hence, the mass is negative for 73 < 1 and for 73 > 2 — again a contradiction. However, 
M > holds in the parameter range 7 G [1,2]. In that case it is not difficult to show that 
the solutions violate the Penrose inequality M > ^J {Ai + A2)/1Qti. Since this inequality is 
related to cosmic censorship, this might indicate that these configurations are not regular 
outside the two gravitational sources. However, since so far no rigorous proof of the Penrose 
inequality for axisymmetric configurations was given, this does not yet exclude the possibility 
that two sub-extremal black holes can be in equilibrium. Instead, one can directly show that 
these solutions always suffer from the presence of singular rings; see Appendix A in (l2| . 
Therefore, this solution branch can be excluded, too. 



VII. SUMMARY 



As a characteristic example for the present discussion about existence or non-existence of 
stationary equilibrium configurations within the theory of general relativity, we have studied 
the question whether two aligned black holes can be in equilibrium. We have shown how this 
question can be reformulated in terms of a boundary value problem for two disconnected 
Killing horizons with specific boundary conditions at the horizons, on the axis of symmetry 
and at infinity. Using the Ernst formulation of the Einstein equations, it was possible to 
apply methods from soliton theory which allow us to study an associated linear problem 
(LP) instead of the non-linear field equations themselves. 
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By integrating this LP along the boundaries of the vacuum region (along the axis, the 
horizons and at infinity) we found an explicit expression for the Ernst potential /"*" on the 
upper part of the symmetry axis A'^ (see Fig. [1]). In particular, it turned out that is a 
quotient of two normalized polynomials in ( (the "axis coordinate"), which depends on a 
number of parameters. In addition, the solution of the LP led to constraints between these 
parameters which we have shown to be equivalent to the "equilibrium condition" (regularity 
condition) a = (a: gravitomagnetic potential) on the axis of symmetry. By analyzing the 
constraints, we found that the Ernst potential /"^ is a quotient of two polynomials of second 
degree in ( and hence all possible equilibrium configurations with two black holes would 
necessarily belong to the double-Kerr-NUT family of solutions. 

Discussing the three possible configurations 

(i) two sub-extremal black holes 

(ii) one sub-extremal and one degenerate black hole 

(iii) two degenerate black holes 

which are characterized by specific restrictions of the parameters ("constraints") of the 
double-Kerr-NUT solutions, we could show that none of them does correspond to a physically 
reasonable black hole solution: All solution families (i), (ii) and (iii) suffer from the presence 
of naked singularities, which manifests in the violation of physical black hole inequalities (the 
positivity of the ADM mass or inequality f ll07p between angular momentum and horizon 
areas for sub-extremal black holes). Hence we can conclude that two-black hole equilibrium 
configurations do not exist. 
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